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SUMMARY 


Research activity in the aerothermodynamics branch at the NASA Ames Research 
Center is reviewed. Advanced concepts and mission studies relating to the next generation 
aerospace transportation systems are summarized and directions for continued research 
identified. Theoretical and computational studies directed at determining flow fields and 
radiative and convective heating loads in real gases are described. Included are Navier- 
Stokes codes for equilibrium and thermochemical nonequilibrium air. Experimental studies 
in the 3.5-ft hypersonic wind tunnel, the ballistic ranges and the electric arc driven shock 
tube are described. Tested configurations include generic hypersonic aerospace plane con- 
figurations, aeroassisted orbital transfer vehicle shapes and Galileo probe models. 


INTRODUCTION 

Current research activity in the Aerothermodynamics Branch at the NASA Ames 
Research Center is outlined and described. Focused activity is directed toward orbiting 
high-altitude, high-drag aerobraking vehicle configurations such as aeroassisted orbital 
transfer vehicles (AOTVs). and high-lift transatmospheric vehicles (TAVs) for the purpose 
of maturing the enabling technology necessary for the design of such vehicles. Research 
efforts are divided into four distinct areas: 1) advanced mission and concept studies, 2) 
computational aerothennodynamic flow-field code development, 3) thermochemical non- 
equilibrium reacting gas models, and 4) code validation experiments. These four research 
areas are considered svnergistically to extend the state of the art in aerothermodynamic 
technology. 

Advanced mission and concept studies identify the thermal-control requirements for 
design-optimized, aeromaneuvering performance for space-based applications and low- 
Earth orbit (LEO, 400 km altitude) sorties involving large, multiple plane-inclination 
changes. Heating analyses for hypersonic, low-density viscous flows including real-gas 
chemistry and surface catalysis are considered. The predicted aerothermodynamic heat- 
ing characteristics are correlated with thermal-control and flight-performance capabilities. 
Mission payload capability for delivery, retrieval, and combined operations is determined 
for round trip sorties extending to polar orbits. The results help to identify technology 
issues required to develop prototype operational systems, and to provide guidance and 
direction for development of computational aerothermodynamic flow codes. 

The design of the next generation aerospace transportation systems will be driven 
by considerations such as fully reusable systems, maximum payload-to-total-weight ratios 
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and high-altitude aeromaneuverability to achieve orbital plane change and cross-range 
capability. These considerations, and others, will be finely tuned to effect designs for 
economically viable vehicles such as rapid-response TAVs and space freighter concepts like 
the AOTV. 

In the past heavy reliance was made using ground-based test facilities such as shock 
tubes, arc jets, and ballistic ranges, in conjunction with engineering design and analysis 
codes, to achieve closure on designs of such vehicles as the Apollo, the Space Shuttle, and 
the Galileo probe. Many of these vehicles were expendable “one-shot” vehicles that use 
ablative heat-shield materials. The space shuttle, a reusable nonaeromaneuvering vehicle, 
employs reusable tiles but must undergo considerable refurbishment between missions. 
To design fully reusable aeromaneuvering vehicles will require design tools significantly 
improved over those used in the past, and design tolerances will need to be much tighter 
to avoid unnecessary weight penalties that can drive the cost beyond economic viability. 

Ground-based test facilities can provide valuable insight for the design and under- 
standing of aerospace vehicles. Such facilities, however, cannot simulate all of the con- 
ditions that will be encountered in planned missions of the future. Shown in figure 1 is 
the flight domain for a variety of aerospace vehicles for typical missions in terms of flight 
Reynolds number and flight Mach number. Superposed are regimes which can be simulated 
by typical ground-based test facilities. The high-altitude, high-speed regime is the most 
difficult to simulate, yet it is this regime that will be encountered by vehicles of the future. 
Flight tests can, and are being planned, such as the aeroassist flight experiment (AFE), to 
provide data in this regime. But these tests are expensive and can provide only a limited 
amount of data. To augment this data and to optimize cost-effective designs, reliance 
must be made on advanced computational techniques. Both ground-based and flight tests 
can provide valuable data to validate these computational methods. The computational 
methods can then, in turn, be used to extrapolate our understanding and analysis in to 
regimes not covered by existing facilities. 

Over the past decade and a half, during the hiatus in hypersonic research, we have seen 
great advances in three key areas that are now being synthesized to develop advanced com- 
putational techniques to a mature technology to enable the design of the next-generation 
aerospace transportation systems. The first of these is the development of inexpensive raw 
computing power. In the early 1970s, scientific computers were of the IBM 360 and CDC 
6600 class. The IBM 370 series and CDC 7600 machines were just becoming available. 
Since that time we have seen the Cray 1 and Cray 2 computers, the CDC Cyber 205 and 
the Fujitsu VP200 and VP400 and Hitachi S 810. Coming along are the ETA 10 and 
the Cray 3 supercomputers. These computers are faster in terms of computing speed and 
larger in terms of high-speed memory by orders of magnitude and, at the same time, cost 
is orders of magnitude less per unit of computing power than a computer of a decade and 
a half ago. 

During the same period a new and powerful tool in aerodynamics was developed; com- 
putational fluid dynamics (CFD). Numerical algorithms have been refined to efficiently 
solve the Reynolds-averaged Navier-Stokes equations for three-dimensional (3-D) steady 
and unsteady compressible flows at transonic and supersonic speeds in an ideal gas. Associ- 
ated with this are sophisticated discretization schemes to describe complex geometries and 
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flow fields. To extend the CFD capability to the hypersonic flight regime experienced by 
aerospace vehicles it is necessary to include “real-gas” effects. Typically this involves the 
description of chemical reactions, thermochemical nonequilibrium and radiative transport 
phenomena. 

At the very high temperatures associated with high-altitude hypervelocity flight there 
is a paucity of experimental data defining collision cross-sections and transition probabil- 
ities necessary to accurately describe the chemical exchange processes involved. To this 
end, recent advances in computational chemistry can be used to fill this critical gap. A 
recent survey describing the capability of computational chemistry and the impact on the 
analysis of orbital transfer vehicles is given by Cooper et al. (ref. l). The combination of 
CFD and real-gas chemistry, along with modern high-speed computers form the synergistic 
basis of computational aerothermodynamics (CAT). Three excellent recent review articles 
on computational aerothermodynamics are given by Anderson (ref. 2), Graves and Hunt 
(ref. 3), and Lewis (ref. 4). 

The particular phenomena addressed by computational aerothermodynamics is the 
prediction of aerothermal loads. These include effects such as aerodynamic forces, con- 
vective and radiative heating rates, gas/surface interactions and surface catalytic effects, 
thermal protection systems (TPSs) with active cooling and plasma layers and their effect on 
communications and power supply. These phenomena are critical to designing vehicles that 
are to be adequately protected from the very high thermal loads that will be encountered 
during hypervelocity flight and that are to be aerodynamically stable and maneuverable. 
Existing ground-based test facilities will be used to perform critical experiments designed 
to provide data for real-gas computer-code validation. 


ADVANCED MISSION AND CONCEPT STUDIES 

In conceptual studies extending over the past two decades (refs. 5-7), the potential 
of aeroassisted technology for enhancing orbital operations and planetary missions has 
been widely recognized. This technique calls for using the aerodynamic forces produced 
by grazing passes through the upper atmosphere to achieve the transition to local orbit 
by deceleration or directional change; earlier methods relied exclusively on propulsion 
power. The propellant saved by eliminating the costly propulsive maneuvers not only 
makes possible missions that are otherwise impractical but also substantially increases 
payload capability. Exploratory studies (ref. 7) have indicated that two classes of AOTVs 
can satisfy a broad range of Earth-centered, space transport missions. 

Aerobraking AOTVs 

One design, useful primarily as a space freighter for transporting large payloads when 
time is not a constraint in the mission requirements, is aerobraking AOTVs. Aerobraking 
vehicles perform its orbit-change maneuvers by aerodynamic drag in the far-outer extent of 
the atmosphere to alleviate surface heat fluxes and pressure forces, therefore minimizing the 
weight penalties for the aeroassist apparatus. Extensive design and mission-performance 
analyses for operations encompassing cislunar space have been given by Menees et al. 
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(refs. 8. 9) and Davies and Park (ref. 10) for an AOTV design in this category. A typical 
near-Earth orbital transfer mission would be between a geosynchronous Earth orbit (GEO, 
35,841 km altitude) and an equatorial LEO. The location of a space station at GEO is 
a future certainty because of its many scientific, commercial, and strategic applications. 
Consequently, low-cost orbital-change maneuvers between GEO and LEO (or Space Shuttle 
orbits) for personnel and material transport will be a future requirement of high frequency. 
The advantages of aerobraking provide substantial savings in propulsion fuel mass for 
this mission. This is illustrated in figure 2 by a typical mission profile for a single-pass, 
AOTV maneuvering between GEO and equatorial LEO. Multiple passes are possible and 
appropriate for unmanned missions for which turnaround time is unimportant, but a single 
pass is appropriate for manned missions. The constraints on the mission are a re-entry 
mass of 12 tons (to conform with other studies) and the use of a liquid rocket engine with 
a specific impulse of 420 sec to make the AOTV reusable and refuelable. The mission 
scenario is as follows: l) the AOTV with initial mass of 30 tons (the approximate capacity 
of Space Shuttle cargo bay) is inserted into equatorial LEO; 2) a propulsive thrust of 2395 
m/sec is required to transfer from the circular LEO to an elliptical orbit with apogee at 
GEO and perigee at LEO; 3) this consumes a propulsion fuel mass of about 45% of the 
initial AOTV mass; 4) a propulsive thrust of 1456 m/sec is required to achieve the circular 
GEO from the elliptical transfer orbit; 5) this maneuver consumes an additional 30% of the 
AOTV mass; 6) a payload of 5.5 tons is picked up for transfer to LEO so that a re-entry 
mass of 12 tons will be obtained subsequently; 7) a retropropulsive thrust of 1490 m/sec 
is required for transfer to an elliptical orbit with perigee in Earth’s atmosphere to take 
advant age of aerodynamic braking; 8) this maneuver expends 30% of the remaining AOTV 
mass in fuel and provides the desired re-entry mass of 12 tons; 9) aerodynamic braking 
occurs with the AOTV, achieving a new elliptical orbit with apogee at LEO; and 10) a 
final small propulsive thrust of 90 m/sec is required for insertion into LEO, which burns 
a fuel mass of only 2.5%. Thus, aerobraking conserves essentially all of the fuel mass that 
would be required for return to LEO using all-propulsive maneuvers, which is about 45% 
of the AOTV mass (the difference between propulsive thrust for steps 2 and 10). 

Several design concepts have been considered for orbital transfer vehicles including a 
conical lifting-brake (refs. 5, 11), a raked elliptic-cone (refs. 12, 13), a raked sphere-cone 
(ref. 14), a lifting body (ref. 15), and a ballute. These five concepts have recently been 
reviewed and critiqued by Park (ref. 16). Historical backgrounds, and the geometrical, 
aerothermal, and operational features of the designs are considered as are the technological 
requirements for the vehicle (namely, navigation, aerodynamic stability and control, after- 
body flow impingement, nonequilibrium radiation, convective heat-transfer rates, mission 
abort and multiple atmospheric passes, transportation and construction, and the pavload- 
to-vehicle weight requirements). These issues are delineated and recent advances that have 
been made in them are summarized. Each of the five design concepts are critiqued and in 
table 1 are rated on these issues. The highest and the lowest ratings are given to the raked 
sphere-cone and ballute designs, respectively. A fairly indepth aeroelastic analysis of the 
bottom-rated ballute concept is given by Park (ref. 17) which identifies critical concerns 
with aerodynamic stability. 
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The highest-rated sphere-cone concept is described in detail by Davies and Park (ref. 
14) and is illustrated in figure 3. To position the thrust line correctly, the two rocket 
engines must exhaust forward through an opening in the aerobrake as seen in figure 3a. 
The thrust line in each of the two engines passes through the approximate center of gravity 
(c.g.) at take-off (fig. 3c) so that, in case one engine fails, there would be no unduly large 
yawing moment. The engines have telescopic extensions that serve two purposes: prevent 
spilling of rocket exhaust gas behind the aerobrake, and increase nozzle area ratio and 
thereby thrust. Two hinged hatches must be provided on the aerobrake for the telescopic 
engines and are indicated in figure '3c. Similar hatches are already incorporated into the 
Space Shuttle for its landing gear. Two sets of fuel tanks exist in this design: the main 
tanks within the heavily armor-protected command/control module, and the auxiliary 
tanks dispersed behind the aerobrake as seen in figure 3b. The fuel in the auxiliary tanks 
will be used first; the fuel in the main tanks (comprising 5-10% of the total fuel capacity) 
will be kept as a reserve. The heavy oxygen tanks are located on the central yaw plane to 
minimize their effect on the trim angle. The tanks for the much lighter liquid hydrogen are 
located symmetrically on the central pitch plane so that possible imbalance in the amounts 
of hydrogen in the two tanks would not greatly affect the trim angle of attack. In case 
one engine fails, oxygen can be shifted from one auxiliary fuel tank to the other to realign 
(in the yaw plane) the c.g. with the thrust line of the functioning engine. The toroidal- 
shaped main oxygen tank inside the command/control module has four compartments. 
Large adjustments in the c.g., if needed during operation of the rocket engines or during 
atmospheric flight, can be made by shifting liquid oxygen among these four compartments. 
Fine adjustments in the c.g. during the atmospheric flight can be performed by gimballing 
the engines. 

On the base side of the aerobrake, the command/control module and the auxiliary 
fuel tanks are covered by a shroud (fig. 3a) made of lightweight heat-shield material. The 
shroud serves three purposes: first, it thermally protects the auxiliary fuel tanks and the 
command/control module from the hot recirculating flow in the base region; secondly, 
it protects the tanks from possible meteoroid bombardment; and thirdly, it is used as a 
radiating surface for the thermal control of the components it encloses. The cargo bay has 
a two-part shield (seen open in fig. 3c) which also serves three purposes. First it thermally 
protects the payload; second, it protects the cargo from meteoroid bombardment; and 
third, in case of failure of one engine, it can serve as a ballast (that is, by opening one 
shield, the c.g. shifts in the yaw plane which brings the c.g. closer to the thrust line of the 
functioning engine). 

High-Lift AO TVs 

The other generic AOTV design is a very high-lift “aeromaneuvering” vehicle. This 
is an essential operational requirement if time-constrained, aeroassisted maneuvers are to 
be accomplished between LEOs involving large, multiple plane-inclination changes. Such 
high-lift vehicles serve as space taxis. They achieve rapid response from one orbital plane to 
another but have the inherent liability of small payload fraction because of low volumetric 
efficiency. 



Results were given by Menees (refs. 7, 9), Davies and Park (refs. 10, 11), and Brown 
(ref. 18) of detailed system-design studies for configurations in the high-lift AOTV cate- 
gory. These were the first studies to address the problem of rarefied-hypersonic flow over 
a lifting surface at incidence with the inclusion of appropriate viscous/inviscid interaction 
phenomena. A subsequent work by Menees (ref. 19) brings together the diverse details of 
the research activities at NASA Ames Research Center for high-lift AOTVs and is an ad- 
junct to this paper. Special emphasis is given to the problems of hypervelocity, low-density, 
viscous effects on flow-field dynamics and thermochemical relaxation effects. 

The use of high-lift vehicles when returning from LEO has many advantages. A 
primary benefit of lift is the ability to achieve large longitudinal range. Lift also provides 
the capability to maneuver within the atmosphere. Maneuver ability enhances mission 
flexibility and increases the choice of landing sites. Lift can also be used to alleviate 
the deceleration loads and aerodynamic heating which accompany atmospheric entries. 
The resultant moderation of the entry environment benefits the vehicle’s occupants and 
permits lower-structural and thermal protection weights. Skipping flight paths can be used 
to maximize both longitudinal and cross ranges (Menees et al. Submitted to AIAA 26th 
Aerospace Sciences meeting, Reno, NV, Jan. 1988) but have associated with them adverse 
deceleration, structural load factors and aerodynamic heating factors when compared with 
gliding flight paths (ref. 20). Even for entry angles of less than 10°, the deceleration loads 
can be an order of magnitude higher for skip flight paths; therefore the crew of a skipping 
vehicle experiences severe physiological stress. To withstand the increased structural loads, 
such major components as the fuselage and wings must be two to three times heavier. 
Heating rates are four to five times higher and require heavier TPSs. 

Air-Breathing TAYs 

The developing need for more economical access to space has spurred renewed interest 
in TAVs in both the United States and Europe. Such a TAV would have responsive, flexible 
operational characteristics approaching those of aircraft and might eventually replace the 
Shuttle by providing short launch notice and turnaround times. The TAV would be able 
to take off and land from ordinary runways and operate in LEOs. A lower-speed variant 
of such a vehicle may be used to transport passengers or high-value cargos to any location 
on the Earth in a small fraction of the flight time of current jet airliners. 

The configurations proposed for TAVs would use primarily air-breathing propulsion 
systems (scramjets) for extended periods of hypervelocity flight to orbital speed. The 
hypervelocity flight within the atmosphere subjects the vehicle to severe aerodynamic heat 
fluxes and total heat loads. Therefore the thermal protection of such vehicles will be major 
design considerations. A TAV using air-breathing propulsion must fly in the denser part 
of the atmosphere to achieve adequate acceleration to reach orbital speed. The elements 
of severe aerothermodynamic environment are therefore coupled with the requirement of 
low aerodynamic drag. To achieve low drag, the vehicle must be slender and must have a 
relatively sharp nose and wing leading edges. Propulsion and aerodynamic performance of 
typical TAVs were approximated and trajectory and angle of attack histories were derived 
which allow meaningful computation of the surface heating rates and heat loads (refs. 21, 
22). Heating rates, equilibrium wall temperatures, and total heat loads at the stagnation 
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point, along a postulated wing leading edge and on the windward surface centerline of the 
vehicle are presented. 

The combination of high heating rates experienced by surfaces with small curvatures 
and the long ascent times results in large total heat loads. Therefore the most severe 
heating occurs during ascent at the stagnation point and wing leading edges. In contrast, 
atmospheric entry occurs at large angles of attack, since high drag is desirable to reduce the 
length of the heating pulse. A comparison of total windward centerline heating loads for 
typical ascent and entry trajectories (ref. 22) are shown in figure 4. Temperatures peak 
around 1300° K during both ascent and entry and radiative cooling should be effective 
over large areas of the vehicle. Ascent-peak stagnation point and wing leading-edge wall 
temperatures are about 3500° K and 2500° K, respectively. Therefore some form of active 
cooling may be required for these regions of the vehicle. The corresponding temperatures 
during entry are 1000° K lower. 


COMPUTATIONAL AEROTHERMODYNAMICS 

Computational fluid dynamics involves the numerical simulation of the equations of 
motion for an ideal gas: these equations are the conservation of mass, momentum and 
energy. In their most general form these equations are the compressible Navier-Stokes 
equations. For turbulent flows the range of length scales is too great to feasibly be resolved 
numerically, and so these equations are time averaged over a scale small with respect to 
the mean motion time scale yet large with respect to the fine-scale turbulent structure. 
These small-scale turbulent transport processes are modelled using eddy viscosity and eddy 
conductivity concepts. Many flow fields have been well simulated for a variety of shapes 
and flow conditions where strong viscous/inviscid interactions and/or flow separation are 
important by solving these equations in a time-like manner until a steady state is asymp- 
totically achieved. When there is no flow reversal and the flow in the streamwise direction 
is supersonic these equations can be simplified by neglecting the streamwise viscous terms. 
The solution to these simplified equations, referred to as the parabolized Navier-Stokes 
equations, can be found by efficient streamwise marching techniques. Further simplification 
can be achieved when viscous/inviscid interactions are weak by decoupling the viscous and 
inviscid dominated regions from one another and simulating the regions separately in an 
iterative manner. Here the inviscid Navier-Stokes equations, termed the Euler equations, 
and are solved in the inviscid region away from body surfaces. Near the body surface the 
viscous-dominated boundary-layer equations are solved. A fourth simplification which can 
be used for strong viscous/inviscid interactions is the viscous shock-layer (VSL) approxi- 
mation. The VSL equations are obtained from the steady-state Navier-Stokes equations by 
retaining terms up to second order in the inverse square root of the Reynolds number. In 
addition, approximations are invoked for the normal pressure gradient and the bow-shock 
location. 

These four equation sets: Reynolds-averaged Navier-Stokes, parabolized Navier- 

Stokes, Euler plus boundary-layer, and VSL, are used in computational aerothermodynam- 
ics to simulate aerothermo loads for a variety of vehicles and flight conditions. Examples 
of each will be described later along with the influence of real-gas effects. 
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Real-gas effects include thermochemical nonequilibrium, where finite-rate processes 
for chemical and energy exchange phenomena occur, and radiative transport. To account 
for chemical reactions, conservation equations for each chemical species must be added 
to the flow-field equation set. There are five flow-field equations; one continuity, three 
momentum, and one energy equation. For dissociating and ionizing air there are typically 
nine species (N 2 , 0 2 , N, O, NO, 0 + , N + , NO + , e~). The inclusion of conservation 
equations for each of these species nearly triples the number of equations to be solved. 
When there are combustion processes or gas/surface interactions or ablation products, 
the number of species increases dramatically. To account for thermal nonequilibrium and 
radiative transport there are additional energy equations to describe the energy exchange 
between the various energy modes (translational, rotational, vibrational, electronic, etc.) 
To further complicate the analysis the range of time scales involved in thermochemical 
processes is many orders of magnitude wider than the mean-flow time scale. This is the 
single most complicating factor in computational aerothermodynamics. A wide variety of 
simplifications are used to alleviate problems associated with widely disparate time scales 
and are discussed briefly next. 

Many flows can be adequately approximated by assuming an equilibrium real gas. 
Here the reaction rates are assumed to all be fast enough that the gas is everywhere in 
local equilibrium and the thermochemical state of the gas can be defined solely by the local 
temperature and pressure. Reactions are allowed to occur but are completely uncoupled 
from the flow-field equations. This is a good approximation for lower altitudes and can 
be used for a major portion of the analysis of such vehicles as transatmospheric vehicles. 
In the other extreme, reactions are sometimes so slow that the gas can be considered 
frozen in a particular chemical state. This phenomenon typically occurs in regions of rapid 
expansion such as in jets or base regions of body shapes, but can sometime be used behind 
compressive shocks as well. 

When finite rate chemical reactions are important they can often be considered to 
be in thermal equilibrium. That is, that the energy modes of the species equilibrate very 
rapidly compared to the chemical rate processes. Even with this simplification the time 
scales vary over an extremely wide range, resulting in a “stiff” behavior of the complete 
equation set and adding to the difficulty in solving the flow-field and species equations in a 
fully coupled manner. In these cases the species equations are often effectively uncoupled 
from the flow-field equations and solved separately in a “loosely” coupled manner, often 
by a different (typically implicit) numerical technique. 

VSL Methods 

Recent applications using VSL techniques are describe in references 18 and 24 - 26. 
Green et al. (refs. 24 . 25) obtained VSL solutions for the forebody of the Titan atmosphere 
entry probe which is planned to make scientific measurements in the organic haze layer of 
Saturn’s largest moon. The flow was assumed to be in local thermochemical equilibrium 
and corresponded to an entry velocity of 12 km/sec. Both radiative and convective heating 
analyses were performed in order to identify the TPS requirements and associated material 
response necessary to protect the probe. Shown in figure 5 is the chemical species distri- 
bution predicted along the stagnation streamline in the forebody shock layer at a time 
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corresponding to the peak heating pulse. The convective and radiative stagnation-point 
heating histories are shown in figure 6. The heating pulse occurs in the hypersonic, low 
Reynolds number region of the trajectory and at the stagnation point the peak flux is 
2.5 MW/m 2 . A stagnation-point convective-heating correlation was developed from the 
results of this study which can be employed in approximate heating analyses. 

Balakrishnan et al. (ref. 26) solved the radiative VSL equations to analyse the flight 
data from several experiments including Fire II, Apollo 4 , and PAET and compared the 
computed results with radiometric data from the flight tests. Calculated convective heating 
rates were found to compare well with experiment, but self-absorbtion and thermochemical 
nonequilibrium effects precluded good agreement with the radiative heating rates. 

In an attempt to account for thermochemical nonequilibrium, Brown (ref. 18) solved 
the VSL equations for dissociating nitrogen for typical flight speeds of 10 km/sec, altitudes 
of 80 km and bodies with nose radii from 1 to 50 cm. A two-species, two-temperature 
nonequilibrium gas model (refs. 27, 28) was assumed and comparisons are made with 
previous theories, experimental data, and results based on the thermodynamic equilibrium 
assumption. Results shown in figure 7, in which heat-transfer rate coefficient is shown as 
a function of rarefaction parameter K r 2 , indicate that the inclusion of vibrational relaxation 
has little effect on the convective heat-transfer rate for a fully catalytic surfaces but can 
contribute to increased heating rates to catalytic surfaces. Preliminary studies (refs. 7, 29) 
indicate that peak heating for a TAV occurs at rarefaction parameter values, K r 2 , between 
8 and 18 for a leading edge radius of 10 cm. If low-catalytic materials are used for thermal 
protection, vibrational nonequilibrium effects will be critical in the structural design. 

The effect of thermochemical nonequilibrium on radiative heating has been addressed 
by Park (refs. 30, 31). The results of these studies indicate that radiative heating loads 
for typical AOTV missions are of the same order as the convective heating loads. 

Euler Methods 

Recent applications using Euler methods are described in references 32-36 (Also 
Palmer, G. E., Submitted to AIAA 26th Aerospace Sciences Meeting, Reno, NV, Jan. 
1988). A total variational diminishing (TVD) numerical scheme for high Mach number 
flows in chemical equilibrium with emphasis on improved convergence is discussed in ref- 
erences 32-34. Here TVD numerical dissipation and real-gas properties are added to one- 
dimensional (1-D) and axisymmetric codes using the Beam and Warming central-difference 
algorithm. Solutions are generated for hemispheres and cylinders at Mach 15 and 20 km al- 
titude. Equilibrium solutions were generated for axisymmetric and two-dimensional (2-D) 
spheres and cylinders at Mach numbers up to 15. The altitude (20 km) and velocities were 
chosen to reflect the flight regime of hypervelocity atmospheric aircraft. The algorithm 
uses a TVD dissipation operator which allows strong shocks to be captured. Real-gas prop- 
erties are provided by the Gordon and McBride (ref. 37) equilibrium chemistry package. 
Particular attention was devoted to achieving quadratic convergence, which no other equi- 
librium code has been capable of duplicating. Density and temperature contours show the 
effect of real-gas chemistry on shock stand-off distance and adiabatic wall temperature. A 
converged solution gives confidence that the method is stable and correct. In the 1-D code 
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(Ref. 32) fully coupled finite-rate chemistry is considered in the simulation of hypersonic 
flow through a strong normal shock. 

Yang (ref. 35) solved the 2-D unsteady Euler equations to simulate the impingement of 
an oblique blast wave on a 2-D AOTV-like configuration. A time-accurate hybrid upwind 
algorithm was used. Shown in figure 8 is the geometry and computational grid used for the 
simulation, and in figure 9 are computed Mach number contours for a Mach 30 impinging 
shock at 15° incidence to the body. These results were computed assuming a frozen gas 
with = 1.1. 

Three-dimensional simulations for equilibrium air are described by Balakrishnan et 
al. (ref. 36). Thermodynamic properties are evaluated at each computational grid point 
using an equilibrium composition method (ref. 37) and the code has been validated through 
detailed comparisons with tabulated data of reference 38. Results were computed for Mach 
numbers of 10, 15 and 20 at entry altitudes of 20 and 50 km for a hemispherical blunt 
body with nose radius corresponding to the Space Shuttle. 

Parabolized Navier-Stokes Methods 

Recent applications using the parabolized Navier-Stokes equations are described in 
references 39-43. Rakich et al. (ref. 39) solved the equations in 3-D for a Space Shuttle 
configuration in ideal air. Balakrishnan (ref. 40) and Prabhu and Tannehill (ref. 41) added 
equilibrium real-air properties. Shown in figure 10 is a comparison of computed temper- 
ature distributions for one cross-sectional plane for both ideal and equilibrium air. The 
real-gas temperatures are significantly lower than the ideal gas temperatures at these flight 
conditions. Prabhu has extended this capability to include chemical nonequilibrium in 2- 
D (ref. 42) and 3-D (ref. 43). The gas dynamic and species concentration equations are 
solved in a coupled manner using a noniterative, implicit, space-marching finite-difference 
method. The conditions for w'ell-posedness of the space-marching method have been de- 
rived from an eigenvalue analysis of the governing equations. Hypersonic laminar flow of 
chemically reacting air over cones and wedges has been simulated and the results vali- 
dated. Shown in figure 11 are computed oxygen and nitric oxide species concentrations in 
the shock layer of a 10° cone for two different altitudes and a speed of 8.1 km/sec. The 
results compare favorably with reacting boundary-layer results (ref. 44). 

Navier-Stokes Methods 

Recent applications using Navier-Stokes methods are described in references 45-50. 
These hypersonic applications are based on two well-established codes. The first, the 
ARC2D/ARC3D codes, have been applied and validated for a variety of complex config- 
urations and flow fields. Equilibrium air properties are included to enhance the codes to 
apply to hypersonic flight. The central-difference Beam- Warming algorithm is enhanced 
by adding TVD modifications so that strong discontinuities can be adequately treated in 
an upwind-like manner. The second, conservative supra-characteristic method (CSCM), is 
an upwind scheme developed for hypersonic applications. Equilibrium air properties are 
included for realistic simulation at flight conditions. Shown in figure 12 are computed con- 
centrations of nitrogen for a Mach 20 flow of equilibrium air over a hemisphere forebody 
(ref. 46). These results compare well with the results of reference 38. 
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Of particular importance in the simulation of complex hypersonic flow fields is the 
ability to adapt the grid to the solution itself. Two codes, ADAPT2D and ADAPT3D 
are especially effective as practical and robust methods to define solution adaptive grids 
(refs. 51-53) These codes have been used for a variety of flows, both steady and unsteady 
and in both 2-D and 3-D. Use of this technique significantly improves the accuracy and 
efficiency of the Navier-Stokes methods. An application of the CSCM method using the 
adaptive grid is shown in figures 13 and 14 for a 2-D hypersonic inlet (ref. 50). The 
solution adapted grid is shown in figure 13 and corresponding computed pressure contours 
on this grid are shown in figure 14. The shock reflection along the lower boundary and 
the interaction of the expansion waves with the reflected shock are the major features of 
this flow. The interaction of the expansion fan with the reflected shock results in a curved 
reflected shock. The predicted pressure recovery behind the reflected shock is critically 
dependent on accurate resolution of these features. Other applications using solution 
adaptive grids with Navier-Stokes methods are described in references 54 and 55, and a 
recent survey of current state of the art in computational aerothermodynamics is given in 
reference 56. 


NONEQUILIBRIUM CHEMISTRY 

The development and validation of a computer program to calculate the flow of nitro- 
gen and air in thermal and chemical nonequilibrium is a first step toward the longer range 
goal of predicting multidimensional flow fields in nonequilibrium air and for hydrogen- 
air combustion processes. This predictive capability is critical to designing the protec- 
tive heat shield for the AOTV where radiation heating loads are expected to be severe, 
and to analyzing scramjet engine performance for the air-breathing propulsion system of 
TAVs. Thermochemical nonequilibrium reacting gas models have been considered by Park 
(refs. 27, 28, 57). A computer program identified as “Shock Tube Radiation Program” 
(STRAP) (ref. 27) has been developed which computes 1-D viscous flows of nitrogen and 
air in a constant-area duct. By improving existing theories and by introducing several 
new innovations, differences between the translational and rotational temperature and the 
vibrational and electron-electronic temperature are considered as are chemical reactions 
whose rates depend on those two temperatures. The computed nonequilibrium thermody- 
namic properties have been analyzed using the Nonequilibrium Air Radiation (NEQAIR) 
(ref. 58) program to determine the spectral radiation characteristics. The resulting radi- 
ation characteristics are compared with available experimental data. Shown in figure 15 
is a comparison of computed and experimental emission intensity spectra for air behind a 
normal shock at a velocity of 10 km/sec and pressure of 0.1 torr. The computations were 
performed using a two-temperature gas model in which the translational and rotational 
temperatures were equal to one another but different from the vibrational and electron 
temperatures which were assumed equal to one another. The point behind the shock 
where the comparison is made is nearly halfway between the onset of the shock and the 
point downstream where the air is once again in thermal equilibrium. Good agreement is 
observed between the experimental and the present theoretical data, thus validating the 
computations and the computer code. 
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Lee (ref. 59) derived the basic governing equations for the low-density, high-enthalpy 
flow regimes expected over the heat shields of proposed AOTVs by combining and extend- 
ing existing theories. The conservation equations are derived from gas kinetic principles 
and a four-component ionized gas consisting of neutral molecules, neutral atoms, singly ion- 
ized ions, and electrons, assuming a continuum flow. The differences among translational- 
rotational, vibrational, and electron temperatures are accounted for as well as chemical 
nonequilibrium and electric-charge separation. Expressions for convective and viscous 
fluxes, transport properties, and the terms representing interactions among various energy 
modes are explicitly given. The expressions for the rate of electron-vibration energy trans- 
fer, which violates the Landau-Teller relation, are derived by solving the system of master 
equations accounting for the multiple-level transitions. Subsequently, a theoretical study 
was made (ref. 60) of the electron-impact vibrational excitation rate processes expected 
in the AOTV flow fields. Semi-empirical quantum-mechanical treatment was adopted to 
predict the vibrational excitation cross sections and the rate coefficients. The obtained 
results for e _ collisions with N 2 showed reasonable agreement with available experimental 
data. The solution to the e“ No vibrational rate equation revealed the unexpectedly slow 
process to equilibrium which is the result of the diffusion characteristic and multiple-level 
transitions at high-temperature regions. A modified Landau-Teller-type rate equation and 
a corresponding relaxation time are suggested which are suitable for numerical calculations 
relevant to AOTV flow fields. 


EXPERIMENTAL PROGRAM 

Experimental facilities used for the purpose of code validation include the Ames 3.5- 
Foot Hypersonic Wind Tunnel which is used for testing generic high-lift configurations 
at Mach numbers of 5. 7. 10 and 14, the pressurized ballistic ranges which can produce 
flight speeds up to 27 km/sec and is used to test AOTVs and planetary probes such as the 
Galileo probe, and the electric arc driven shock tube which can produce nonequilibrium 
flows at speeds up to 13 km/sec. Flight experiments for both the drag brake and the 
high-lift configurations are planned for early in the next decade. 

Shock Tube 

The 24-in. EAST Facility w’as operated with pressures between 7 and 30 gm of mer- 
cury to simulate flight altitudes between 68 to 81 km. Test firings show that: l) shock 
velocities in excess of 13 km/sec are achieved at densities equivalent to altitudes of 80 km., 
2) the hot driver gas emits radiation in the same pattern as observed in previous higher 
density, slower shock speed tests, 3) the driver gas does not radiate at unanticipated spec- 
tral frequencies, and 4) there is no significant radiation spectra from common impurities 
such as sodium and iron. This facility will be used to perform basic experiments on ther- 
mochemical nonequilibrium with particular emphasis on vibrational energy excitation and 
nonequilibrium. A description of proposed tests and instrumentation is given by Sharma 
and Park (ref. 61). 
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Ballistic Ranges 

Two ballistic range facilities at Ames Research Center support hypersonic research 
(ref. 62), the Hypersonic Free Flight Aerodynamic Facility (HFFAF) and the Pressurized 
Ballistic Range (PBR). The two facilities compliment each other; each has specific advan- 
tages for certain types of tests. Together, they allow testing of a variety of models over a 
wide range of test conditions. Each facility and associated experiments is described briefly. 

Hypersonic Free Flight Aerodynamic Facility 

A photograph and schematic of this facility are shown in figure 16. The test section is 
23-m long and has 16 orthogonal spark shadowgraph stations evenly spaced (1.52 m) over 
its length. Kerr-cell shutters are used to produce a sharp model and flow-field image on 
the film. Four deformable-piston light-gas guns, having bore diameters of 0.71, 1.27, 2.54, 
and 3.81 cm, are available for launching the model into free flight. Each of these guns can 
operate to muzzle velocities of about 9 km/sec. A shock tube to provide a countercurrent 
flow capability is available, but has not been used for recent tests. 

Tests in this facility can be conducted from 1 atm to as low as about 20 pm Hg. And 
tests can be conducted in nontoxic gases other than air, such as CO 2 , H 2 , He, Kr, and Xe. 

Pressurized Ballistic Range 

The PBR consists of a large tank that can be pressurized or evacuated and a test 
section 62-m long. It has 24 orthogonal spark shadowgraph stations irregularly spaced 
over its length. The station spacing ranges from 2.1 m to 4.2 m. All of the optics are 
internal to the tank, which imposes a limit on the maximum velocity of the model. Many 
tests in this facility are conducted with powder gas guns rather than the more sophisticated 
deformable-piston light-gas guns. The advantages of this facility over HFFAF are threefold: 
long model trajectory, pressures up to about 6 atm, and highly detailed shadowgraphs. 

Galileo Tests 

The ballistic ranges at Ames have supported all of the U.S. probe missions to other 
planets. These include the 1976 Viking Mission to Mars and the 1978 Pioneer Venus 
Mission. Currently, tests are being conducted on the Galileo probe. The Galileo spacecraft 
will be launched in the near future and will arrive at Jupiter 2 yr later. The spacecraft 
consists of two major components, an orbiter which is to orbit Jupiter numerous times, 
concentrating on close encounters with many of Jupiter’s moons, and a probe to enter (and 
descend through) the atmosphere of Jupiter. The probe will make in situ measurements 
as it descends through the atmosphere prior to its eventual destruction due to extreme 
external pressures. 

Although the probe aerodynamics were needed initially for design purposes, more 
accurate aerodynamics are needed in support of the Atmosphere Structure Experiment 
carried on board the probe. The Atmosphere Structure Experiment is designed to deter- 
mine the state properties (i.e., density, pressure, temperature) of an unknown planetary 
atmosphere as functions of altitude from measurements made during the entry and descent 
of a probe. The experiment consists of a three-axis accelerometer, plus pressure and tem- 
perature sensors. During the high-speed portion of the trajectory, from an entry velocity 
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above 47 km/sec to sonic speed, direct measurements are impractical and accelerometers 
are used to determine the state properties. This requires the precise knowledge of the 
probe aerodynamics, in particular the vehicle drag and lift coefficients as functions of 
Mach number and Reynolds number. The aerodynamic characteristics plus the measured 
decelerations allow the probe attitude to be determined and the atmospheric density to 
be deduced. Integration of the density gives the pressure, and the temperature is deduced 
from the equation of state (given the molecular weight, which is measured by another on- 
board experiment). The ballistic range facilities are well suited for providing the accurate 
aerodynamic data over a wide range of conditions. 

The Galileo probe is a blunt, 45° cone. During the high-speed part of the entry, 
severe ablation takes place. As much as 40% of the vehicle mass at entry is expected to 
be ablated away, mostly in the nose and conical regions, and even the maximum diameter 
will be significantly decreased. Hence, tests must be conducted not only of the entry 
configuration but of hypothetical ablated configurations as well. 

A typical shadowgraph of a Galileo model in flight obtained in HFFAF is shown in 
figure 17. The screw on the model base is for attachment to its sabot prior to launch. 
The vertical wires are plumb bobs on both sides of the facility for reference purposes and 
the irregular markings are imperfections in the facility windows caused by past impacts of 
various sorts. 

Low Reynolds Number Tests 

Ballistic range tests are currently being conducted to precisely define the drag char- 
acteristics of the Galileo probe at Reynolds numbers, based on model diameter, of about 
500 to 250. The importance of obtaining drag data at these low Reynolds numbers is 
due to the fact that the drag coefficient is expected to increase markedly as the slip-flow 
and free-molecule-flow regimes are approached. This dramatic increase in drag coefficient 
occurs below a Reynolds number of about 1000. 

The importance of obtaining drag data at various Reynolds numbers is shown in 
figure 18. Shown are Pioneer Venus data down to a Reynolds number of about 250. The 
drag coefficient increases continuously below a Reynolds number of 1 million, but the 
increase becomes most dramatic below 1000. 

AOTV Tests 

Ballistic range tests have been conducted for two AOTV configurations; a symmetric 
configuration and raked elliptic cone. These tests were conducted to: 1) Provide ex- 
perimental aerodynamic data and good flow-field definition against which computational 
aerodynamicists could validate their computer codes. 2) Define bow shock-wave shape 
and shock standoff distance, 3) Investigate flow impingement on afterbody, 4) Compare 
the aerodynamics of several configurations, 5) Investigate how minor changes in corner 
geometry affect the flow field, 6) Determine the trim angle of attack of a trimmed vehicle. 

A shadowgraph for the symmetric configuration in figure 19 and for the raked elliptic 
cone in figure 20. These flow visualizations plus drag data from these tests are used for 
computer code validation. 
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3.5-Ft Hypersonic Wind Tunnel 

The current test program in the 3.5-ft hypersonic wind tunnel is focused on a generic 
all-body (elliptic cross section, delta planform) hypersonic aircraft model. This config- 
uration is representative of airbreathing TAVs now being considered. Flow-visualization 
(shadowgraphs and surface oil-flow patterns), surface-pressure, surface heat-transfer, and 
flow'-field surveys using probes and nonintrusive lasers will be obtained for this model both 
with and without control surfaces. Angles of incidence to the free stream will be varied 
between plus and minus 15° and the free-stream Reynolds number varied between 1.5xl0 6 
and 25xl0 6 . Shown in figure 21 is a model with all control surfaces attached. The model 
has a half angle of 75° and is 1-m long. There are two nose configurations: a sharp nose 
and a blunted nose. The simple geometrical configuration is easy to define for computer 
flow codes and can be readily used to evaluate a wide variety of hypersonic flow codes. 

Aeroassist Flight Experiment (AFE) 

A forthcoming NASA flight experiment called the AFE is planned for early in the 
next decade. Ames Research Center will participate in this experiment by assimilating a 
large base of radiometric data for high-altitude, high-velocity, thermochemically nonequi- 
librated flow conditions. The AFE will be carried to orbit by the space shuttle and then 
deployed for the atmospheric data pass. Accelerated by an 18,000-lb-thrust solid-rocket 
motor, the vehicle will enter the atmosphere at nearly 10 km/sec, and then experience 
approximately 500 sec of aerodynamic deceleration during which a variety of flight data, 
including radiative and convective heating rates, will be gathered. The vehicle will exit the 
atmosphere at orbital speed to be recovered by the shuttle orbiter for return to earth for 
post-flight evaluation. As a preliminary to the design of a radiometer for this experiment, 
an approximate method for predicting both equilibrium and nonequilibrium radiative sur- 
face fluxes has been developed (ref. 63). Spectral results for one trajectory state, a velocity 
of 10 km/sec at an altitude of 85 km, are shown in figure 22, where the spectral surface 
flux at a distance of 20.9 cm behind the shock front is plotted as a function of wavelength 
in the spectral region from 0.2 to 2.0 pm. An inspection of the figure reveals that the 
spectrum appears to be composed of a background continuum with a color temperature 
in the range of 7.000° to 8,000° K (based on a flux maximum in the vicinity of 0.4 pm) 
on which is superimposed a complex structure of molecular bands and broadened atomic 
lines. The radiation calculation included nine species (O 2 , N 2 , NO, O, N, N + , 0 + , N T 2 ^, 
e“); some of more apparent band-heads and lines from these species are identified in the 
figure. These results, and others like them, are used to develop the instrument parameters 
for the three different types of radiometers proposed for the experiment. 


CONCLUDING REMARKS 

Aerothermodynamic research underway at NASA Ames has been described. Four re- 
search areas are considered synergistically in an effort to mature the enabling technology 
necessary for the design of the next generation aerospace transporation systems. These 
four areas include 1) advanced mission and concept studies, 2) computational aerother- 
modynamic flow-field code development, 3) thermochemical nonequilibrium reacting gas 
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models, and 4) code-validation experiments. Research in these areas continues in an effort 
to improve the accuracy and efficiency of predictive methods and our understanding of 
hvpervelocity flows and its effect on aerospace vehicles. 
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Figure 1.- Flight-domain simulation capability. 
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Figure 2.- Mission schematic for AOTV in coplanar maneuver between LEO and GEO. 
(Masses computed for ip = 420 sec.) 
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Figure 4.- Total heat loads on generic TAV windward centerline. 
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Figure 5.- Heating histories at the stagnation point (S/R n = 0). 



Figure 6.- Shock-layer structure during peak heating along stagnation streamline - chem- 
ical species profiles. 



Figure 7.- Convective heat-transfer rate correlation. 
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Geometry and grid for two-dimensional AOTV flow-field simulation. 
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Figure 11.- Computed species mass-fraction profiles in the shock layer of a 10° cone. 
Voo = 8.1 km/sec. a) Oxygen mass fraction; b) Nitric oxide mass fraction. 
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Figure 12.- Comparison of N 2 mole-fraction contours over hemisphere nose, 
a) CAG12 solution; b) Solution from reference 38. 





5 . 0 . 


32 




ORIGINAL PAGE IS 
OF POOR QUALITY 

Poo = 0.1 torr 



Figure 15.- Comparison between calculated and measured spectra behind a plane normal 
shock wave. 




Figure 16.- Hypersonic Free Flight Aerodynamic Facility. 
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Figure 17.- Shadowgraph of Galileo Probe. = 14, Re<* = 100,000. 


Figure 18.- 



Effect of Reynolds number on drag-character characteristics of Pioneer Venus. 
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Figure 19.- Shadowgraph of symmetric AOTV. 



Figure 20.- Shadowgraph of AFE model. 
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Figure 21.- Generic all-body hypersonic aerospace plane model, 
with elliptic cross section. 
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Figure 22.- Nonequilibrium-flow spectral flux: Shock-layer Depth = 20.9 cm. 
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